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ABSTRACT 



We describe a new approach to fitting the UV-to-optical spectra of B stars to model atmospheres and present 
initial results. Using a sample of lightly reddened stars, we demonstrate that the Kurucz model atmospheres can 
produce excellent fits to either combined low dispersion IUE and optical photometry or HST FOS spectrophotom- 
etry, as long as the following conditions are fulfilled: 



1 . an extended grid of Kurucz models is employed, 



2. the IUE NEWSIPS data are placed on the FOS absolute flux system using the Massa & Fitzpatrick (1999) 

transformation, and 
> ■ 

3. all of the model parameters and the effects of interstellar extinction are solved for simultaneously. 

When these steps are taken, the temperatures, gravities, abundances and microturbulence velocities of lightly 
reddened BO-AO V stars are determined to high precision. We also demonstrate that the same procedure can be 
used to fit the energy distributions of stars which are reddened by any UV extinction curve which can be expressed 
by the Fitzpatrick & Massa (1990) parameterization scheme. 

We present an initial set of results and verify our approach through comparisons with angular diameter mea- 
rS^ surements and the parameters derived for an eclipsing B star binary. We demonstrate that the metallicity derived 

Q-f from the ATLAS 9 fits to main sequence B stars is essentially the Fe abundance. We find that a near zero mi- 

croturbulence velocity provides the best-fit to all but the hottest or most luminous stars (where it may become a 
surrogate for atmospheric expansion), and that the use of white dwarfs to calibrate UV spectrophotometry is valid. 

Subject headings: stars: atmospheres — stars: early-type — stars: fundamental parameters — stars: abundances 
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1. INTRODUCTION surface gravity logg). Second, the shape of the B star SED's 

>J> T , , ,.„ , r ., , ■ . , j. .., in the UV is also influenced by opacity due to Fe group ele- 

X] ■ The goal of modeling the observed spectral energy distribu- ._ _, „ _ .... , , , , , „ . 

t , ■ .• /or^\ c * • Z • x.te j i * u • • * ments (Cr, Mn, Fe, Co, Ni), and thus should allow constraints 

*H tion (SED) of a star is straightforward, namely, to obtain infor- \ ' _ , 

., , . , ... c . j. ., . , , to be placed on Fe group abundances — without the model- 

. . . i mation on the physical properties of a star from the model at- . .",..,,,. b *1, _ , , . , , ... 

, l'lc, j. l .- u » Aft 4.u u t c^- mg of individual line profiles. Such results, particularly if at- 
mosphere which fits the observations best. After the best-fitting ? . , .,, *,. , ,,, , ,- , 
j i , • c j •* • j ^u i^u 4 u- u a tamable with low -resolution data, would be extremely useful 
model is found, it is assumed that the parameters which de- , , . „ „. . , , • , ■ , 
c ., , , , cc j. . „ , tor studying stellar metallicity patterns when combined with 
fine the model (effective temperature, surface gravity, elemen- ,, , , • • r ^,^ , , n 
. , , ■ . n e ... e „ ... ,,, . •. , e ground-based determinations of CNO abundances (e.g., Smartt 
tal abundances, etc.) faithfully represent those of the star itself. = _, ,, ,„„„ ^ r ,^^o ^ , , ,^^on 
rr. , , ',. /,,. , *• • vaa a & Rolleston 1997, Dufton 1998, Gummersbach et al. 1998) 
The degree to which this key assumption is valid depends upon . , , , _ , _ „.,..,/ 
, ,, ., . c ., c , , ■ ■ c a 4. a- c since only the coolest of the B stars have sufficient lines in their 
both the uniqueness of the fit (requiring a firm understanding of . , J ,, _ , , 
„ , ... , , -v. • t c 4i. u • i optical spectra to allow accurate Fe group abundances to be de- 
all observational errors) and the appropriateness of the physical • , ,- , , 

rived from the szround 

assumptions used to construct the models (which often can be . . , 6 . . „ , _, . , 

■j- j, •„ • e ,. s „., r.. A . Accessing the information content of the B star UV/optical 

verified by ancillary information). It the fit provides strong con- „,™, . , , , ,. , , 

..... j, . j ., „ a i • -a a SED s requires a set of model atmospheres which reproduce 

straints on the model parameters and the model is considered , , • , , % r , „ 

,. , ., .!.•.< t - £„■ „ • . c i the physical structure and emergent flux of the B star atmo- 
vahd, then this continuum htting process provides a powerful . * J __ . _ b , , 
. , c . , ■ . ,, 4- spheres. The main sequence B stars are the hottest, most mas- 
tool for studying stellar properties. • , , • , • 

r™ . ■ i. .. c a v ,u TTiTi 4- i sive, and most luminous stars whose atmospheric structures are 

There are two maior motivations for modeling the UV/optical , ,. , „ ,,,-,,-• 

(, rn , £ ., • n . ff , . %„. mnro c believed to be well-represented by the simplifying assumptions 

SED s of the main sequence B stars. First, since 70-to-100% of ., , , , *\ .... r /T1 ™i^ , ,, , 

the energy of B stars (with 10000 < T eff OOOOO K) is emit- of local thermodynamic equilibrium (LTE), plane parallel ge- 

11 ometry, and hydrostatic equilibrium. Thus, although it is clear 

ted in the accessible spectral region between 1216 A and 1 ^ a m non . LTE analysis may be needed for detailed model . 

M m (see Fig. 1 of Bless & Percival 1997), even low-resolution ing of ific Hne files> LTE modd atmosphere ca i cu i ations 

UV/optical observations can place strong constraints on stellar m expected to re p roduce the atmospheric structures and gross 

models and should provide robust determinations of the stellar emergent flux distributions of these stars (e.g., Anderson 1985, 
properties which influence the SED (particularly T e ff and the 



1989; Anderson & Grigsby 1991; Grigsby, Morrison, & An- 
derson 1992; Gummersbach et al. 1998). As a result, the line- 
blanketed model atmospheres computed by R.L. Kurucz, which 
incorporate the simplifying assumptions listed above, have be- 
come a standard tool in the study of B stars. These ATLAS - 
series models have been used in a wide variety of investigations, 
including the determination of fundamental stellar parameters 
(e.g., Code et al. 1976); population synthesis (e.g., Leitherer 
et al. 1996); stellar abundances within the Galaxy (e.g., Geis 
& Lambert 1992; Cuhna & Lambert 1994; Smartt & Rolleston 
1997; Gummersbach et al. 1998) and among the nearby Galax- 
ies (e.g., Rolleston et al. 1996 and references therein); the phys- 
ical structure of atmospheres of the pulsating (3 Cephei stars 
(e.g., Cugier et al. 1996); the derivation of interstellar extinc- 
tion curves (e.g., Aannestad 1995); the determinations of the 
distances to binary Cepheids (e.g., Bohm-Vitense 1985); and 
many others. 

Despite the long history of "continuum fitting" the 
UV/optical SED's of B stars (see, e.g., Code et al. 1976, Under- 
bill et al. 1979; Remie & Lamers 1982; and Malagnini, Farag- 
giana, & Morossi 1983; among many others), several key is- 
sues regarding the applicability of the ATLAS models to the B 
stars have not yet been addressed quantitatively. These include 
the range of stellar properties over which the physical assump- 
tions and simplifications are valid (or, at least, not debilitating) 
and whether — within that range — the models truly repro- 
duce the observations. Most previous applications of the con- 
tinuum fitting technique have relied on "by eye" fitting, have 
attempted to deduce only T e ff, and have utilized crude treat- 
ments of interstellar extinction. Thus it is not clear whether the 
often large discrepancies seen in the UV region, in particular, 
are attributable to failures in the models or to shortcomings in 
the fitting process. 

Because of their ubiquity, great diagnostic potential, and rel- 
ative ease of calculation, the ATLAS models are a tremen- 
dous resource for stellar astronomy and a quantitative assess- 
ment of their range of applicability is imperative. Recent sig- 
nificant advances in the requisite theoretical and observational 
material make this an ideal time to undertake such an inves- 
tigation. These advances include the release of the ATLAS 9 
models with their greatly improved opacities (Kurucz 1991), 
the re-processing of the entire International Ultraviolet Ex- 
plorer (IUE) low-resolution UV spectrophotometry archives, 
the derivation of a precise UV calibration by the HST Faint 
Object Spectrograph (FOS) team (Bohlin 1996, and see Bless 
& Percival 1997), the derivation of a rigorous transformation 
between the IUE and FOS calibration systems (Massa & Fitz- 
patrick 1999, hereafter MF99), and an improved understand- 
ing of the properties of UV/optical interstellar extinction (Fitz- 
patrick & Massa 1990, hereafter FM; Fitzpatrick 1999, here- 
after F99). 

This paper is the first in a series whose goals are to quan- 
tify the degree of agreement between the ATLAS model at- 
mospheres and UV/optical spectrophotometry of B stars and 
to apply the models to determine the physical properties of the 
stars and the wavelength dependence of interstellar extinction. 
Here we discuss the methodology to be employed by our study 
and present some initial results. In §|[ we provide an overview 
of modeling observed energy distributions and describe the es- 
sential ingredients of the fits: a complete set of model atmo- 
sphere calculations, and a quantitative characterization of the 
effects of interstellar extinction. The section concludes with a 
description of the mathematical procedure used to fit the models 



to the observations and a brief discussion of the determination 
of metallicity using the ATLAS 9 models. In §|| we describe 
the spectrophotometric observations we will be modeling and 
§|] provides a sample of our results, illustrating the quality of 
the fits, the agreement between the derived parameters and fun- 
damental measurements, and a comparison with previous, in- 
dependent abundance measurements. The potential for future 
study is also discussed. Finally, §J5| summarizes our findings 
and outline a few remaining technical issues. 

2. MODELING THE OBSERVED STELLAR FLUXES 

In this section, we provide an overview of the general prob- 
lem of fitting stellar energy distributions, describe the con- 
stituents of the model used to fit the observations, and discuss 
how the fitting is accomplished. 

2.1. Overview of the problem 

The energy distribution of a star as received at the earth, f\$, 
depends on the surface flux of the star and on the attenuating ef- 
fects of distance and interstellar extinction. This observed flux 
at wavelength A can be expressed as: 



h =F<!ni).,x (^) xio- 1 ^' 



(1) 



= F({a}) x x f ^ ) x 10-°- 4£(B - WA - V)+S(l/)] (2) 

where F\ is the emergent flux at the surface of the star (as- 
sumed to be single), {a} is the set of intrinsic physical proper- 
ties which determine the emergent flux from the star, R is the 
stellar radius, d is the distance to the star and A(A) is the total 
extinction along the line of sight. The factor (R/d) 2 is essen- 
tially equivalent to the square of the angular radius of the star 
and is often written as (8/2) 2 , where is the stellar angular 
diameter. 

In eq. 0, the extinction terms are rewritten as quantities 
normalized by E(B-V), i.e., the normalized extinction curve 
k(\-V) = [A(X)-A(V)]/E(B-V), and the ratio of selective to 
total extinction in the V band R(V) =A(V)/E(B-V). 

2.2. The models 

We represent the stellar surface fluxes with R.L. Kurucz's 
ATLAS 9 line-blanketed, LTE, plane-parallel, hydrostatic 
model atmospheres. These models are functions of the 4 pa- 
rameters {a} = {T e ff,logg,[m/H],v t } - which are the effec- 
tive temperature, the log of the surface gravity, the logarith- 
mic metal abundance relative to the Sun, and the magnitude 
of the microturbulent velocity field, respectively. The standard 
ATLAS 9 calculations provide emergent fluxes at 1221 square- 
binned wavelength points from 90.9 A to 160 /i. The bin widths 
are typically 10 A in the UV (300 A < A < 3000 A) and 20 A 
in the optical (3000 A< A < 10000 A). 

Because the fitting procedure must be given access to the en- 
tire region of parameter space which is physically reasonable, it 
was necessary extend the standard Kurucz grid, which contains 
only v, = 2 km s" 1 models. This was accomplished using the 
ATLAS 9 programs and input Opacity Distribution Functions 
(ODFs) to compute additional model atmospheres and emer- 
gent fluxes. We expanded the grid to include v, =0,4,8 km s , 
[m/H] =-1.5, -1.0, -0.5, 0.0, +0.5 (i.e., from 1/30 to 3x so- 
lar), 36 values of T eff between 9000 and 50000 K, and typically 



~ 10 values of logg for each T e ff extending from the Edding- 
ton limit to logg = 5.0. Our final grid consists of emergent flux 
distributions for about 7000 models. 

In order to calculate a model at a set of parameters 
{T*ff,logg*, [m/H] ,v*} which do not exactly coincide with 
a grid point, an interpolation scheme in all 4 parameters is re- 
quired. We adopt the following approach: First, at the three 
grid values of v t closest to v* and the three grid values closest 
to [m/H]*, we produce models with T*,-,- and logg* through bi- 
linear interpolation of \ogF\ against logT e ff and logg. This 
yields 9 models. Next, we quadratically interpolate \ogF\ 
against [m/H] at each of the v,, producing 3 models of the 
desired T*tp logg*, and [m/H]*. Finally, we quadratically 
interpolate log Fa against v, to obtain the final model. Tests 
of this procedure against models computed "from scratch" at 
{T*f t, logg* , [m/H]* , v* } indicate that the errors introduced are 
negligible. 

In addition to fluxes, we produced synthetic Johnson UBV 
and Stromgren uvby photometry for each model, using the pro- 
grams of Buser & Kurucz (1978) and Relyea & Kurucz (1978). 
The value of the synthetic V magnitude V syn — which is used in 
the current fitting procedure — is determined by the following: 



V syn = -21.07-2.5 x log / f x Sx(V)d\ 



(3) 



where f\ is a stellar SED (either observed or theoretical) and 
S\ (V) is the normalized filter sensitivity function for the V fil- 
ter from Azusienis & Straizys (1969). The integration is per- 
formed over the V filter's range of response (4750 A to 7400 
A). The normalization zeropoint, -21 .07, was determined from 
synthetic photometry performed on FOS spectrophotometry of 
seven calibration stars (Colina & Bohlin 1994; Bohlin, Col- 
ina, & Finley 1995) and on ground-based spectrophotometry of 
HD 172167 (Vega; Hayes 1985). For these eight stars, equa- 
tion 3 yields a standard deviation of 0.01 mag in the difference 
V - Vgyn, consistent with expected observational error. 

2.3. Interstellar extinction 

Correcting stellar SED's for the effects of interstellar extinc- 
tion is problematical because the wavelength dependence of ex- 
tinction is known to be highly variable spatially. While aver- 
age extinction curves have been defined, there may be few if 
any sightlines for which some globally-defined mean curve is 
strictly appropriate. Furthermore, the error introduced by using 
an inappropriate curve increases in proportion to the amount 
of extinction (see Massa 1987). Consequently, we adopt one of 
two approaches for modeling extinction, depending on the color 
excess of the object. For lightly reddened stars (E(B-V) < 0.1) 
we adopt a universal mean curve and solve only for E(B-V). 
For the more reddened stars, we investigate the feasibility of 
deriving both E(B - V) and the shape of the extinction curve 
from the data themselves. In both cases we utilize an analytical 
representation for the extinction wavelength dependence based 
on the work of FM. 

For the universal mean curve, we use the results of F99. 
These are based on the general UV curve parameterization of 
FM and the discovery by Cardelli et al. (1989) that the shape 
of UV extinction curv es correlate with the extinction parameter 
R(V), defined in §2.1 F99 presents a mean curve for the case of 
R(V) = 3.1 which (when applied to a stellar energy distribution 
and convolved with the appropriate filter responses) reproduces 
the mean extinction ratios found in the optical and IR from the 



Johnson and Stromgren photometric systems. The curve also 
joins smoothly with the UV region (at A < 2700A) via cubic 
spline interpolation. The UV portion of the curve is expressed 
as a set of three smooth functions and specified by 6 parameters, 
all adjusted to their nominal values for the case R(V) = 3.1. 

For the more heavily reddened stars, we allow for additional 
free parameters (i.e, in addition to E(B-V)) so that a "cus- 
tomized" extinction curve can be derived by solving for some 
or all of a modified version of the FM curve parameters. FM 
showed that any normalized UV extinction curve can be rep- 
resented by 6 parameters which describe the general slope and 
curvature of the curve and the position, width, and strength of 
the 2175 A bump. We modify this approach slightly, in that the 
linear terms are combined (FM show that these are functionally 
related) and we use a cubic spline to bridge the gap between 
the UV at 2700 A and the normalized optical curve near 4100 
A (see F99). Consequently, the most general UV/optical ex- 
tinction curve contains 5 free parameters. The feasibility of this 
approach is demonstrated in §[|. 

2.4. Fitting the spectra - qualitative considerations 

The ability to extract the physical parameters of a model ac- 
curately from the observed SED depends on whether changes 
in those parameters produce detectable and unique spec- 
tral signatures. In Figure |J we illustrate the signatures of 
the five principal parameters considered in the fitting pro- 
cess for a typical B star model flux distribution. The top 
panel shows the UV/blue-optical SED for an ATLAS 9 model 
with parameters as listed in the figure. The five lower pan- 
els show the effects on the spectrum caused by changes 
in the physical parameters {AT e ff, Alogg, A[m/H], Av,} = 
{+100K,+0.1dex,+0.1,dex,+lkms" 1 }andbyachangeA£'(B- 
V) = +0.005 mag in the amount of interstellar extinction 
present. 

Three aspects of the figure are important: 

1 . The spectral signatures of the four stellar parameters are 
functionally independent, i.e., it is not possible to con- 
struct any one effect from a linear combination of the 
others. Consequently, it should be possible to derive un- 
ambiguous and uncorrelated estimates for each of the pa- 
rameters from the data. 

2. The signature of interstellar extinction is quite distinct 
from those of the model parameters and, in principal, the 
value of E(B - V) should be precisely determined by the 
fitting procedure. 

3. The RMS magnitudes of each of the effects in the UV re- 
gion are of order 0.01 dex. Since we will be fitting more 
than 100 data points in the UV with typical RMS errors 
of 3% per point, uncertainties in the physical parame- 
ters corresponding to the differences shown in the figure 
should be well within our grasp. 

Note that the strength and shape of each of the spectral sig- 
natures shown in Figure [I] — and their degree of independence 
from each other — do vary as a function of the model param- 
eters, particularly T e ff. Therefore, the precision achievable in 
estimating the parameters will differ across the wide range of 
T e ff values appropriate for the B stars. 

It might seem surprising that v, and [m/H] can each be de- 
termined independently of the other, since both directly affect 



the absorption lines. This can be understood in principle, how- 
ever, from considering the curve-of-growth, which describes 
the strength of an absorption feature as a function of the various 
relevant physical properties of a system. In particular, varia- 
tions in the number of available absorbers (e.g., [m/H]) mainly 
affect weak, unsaturated absorption lines on the "linear part" of 
the curve-of-growth, while variations in line broadening (e.g., 
v t ) mainly affect the stronger, saturated lines on the "flat part" 
of the curve-of-growth. Thus, changes in [m/H] and v, have 
their greatest effects on distinct populations of lines (although 
there is a large "middle ground" in which both parameters are 
important). This can be seen in Figure [l] where the spectral fea- 
tures near 2000 A respond dramatically to a change in v, and 
much less so to a change in [m/H]. The inverse can be seen for 
features in the 1300-1500 A region. 

The ability to distinguish microturbulence effects from abun- 
dance effects appears to breakdown for cases where v, is large. 
In §4 we show several stars for which the combination of a high 
microturbulence and a low metallicity reproduce the observed 
UV opacity very well — although their individual values are 
unlikely to represent the corresponding stellar properties. As 
discussed in §4, this probably indicates that microturbulence is 
a physically incorrect description of the line broadening mech- 
anism in these stars. 

2.5. Fitting the spectra - quantitative considerations 
In practice, we fit equation (g) recast into the following form 



log /a© = log 



-) xl0-° 4,!Vi 
d 



+logf(7; // ,logg, [m/H], v,) A 
-0AE(B-V)k(\-V) 



(4) 



We refer to the leading term as the "attenuation factor" because 
it is independent of wavelength. It contains the distance, stel- 
lar size and total extinction information. The second term is the 
stellar model which has a non-linear wavelength dependence on 
the 4 model parameters. The final term expresses the magnitude 
and wavelength dependence of interstellar extinction and it can 
represented as either a simple linear function (in the case of pre- 
defined extinction curve) or a more complex function involving 
both linear and non-linear terms (see §2.3). In general, there 



can be anywhere between 6 and 1 1 free parameters, depending 
on how much flexibility is incorporated into the description of 
reddening. 

Notice that, since we do not utilize infrared photometry in 
our program, we cannot determine the value of R(V) from the 
fitting procedure. Instead, we adopt a typical value (usually 
R{V) = 3.1) and then incorporate the expected uncertainty into 
the error estimate of log(R/d) 2 (an uncertainty in R(V) does not 
affect the stellar parameters or k(X-V)). 

The best-fitting model is obtained by solving equation (Q) 
for all of the unknown parameters simultaneously. We use an 
iterative, weighted non-linear least squares technique (the Mar- 
quardt method) to adjust the parameters and minimize the \ 2 of 
the fit. The parameter values used to initialize the fitting proce- 
dure are determined from the spectral type and optical photom- 
etry of the star. The only constraint we impose on the domain 
of the model parameters is that v, and E(B-V) be > 0. Note 
that E(B-V) is simply a scaling factor for k{X-V), and we do 
not utilize (B — V) photometry in the fits at this time. 



The degree of accuracy to which any single parameter can be 
determined depends on the sensitivity of the observed flux to 
that parameter and on the degree of covariance with the other 
parameters (e.g., see Fig. |j). Our error analysis incorporates 
the full interdependence of the various parameters and includes 
the effects of random noise in the data. All uncertainties quoted 
in this paper are \-a internal fitting errors, in the sense that a 
change in a particular parameter by ±l-<7 results in a change 
in the total x 2 by +1 a ft er all the other parameters have been 
reoptimized (see Bevington 1969, Chapter 11). 

2.6. Abundances from ATLAS 9 models 

The definition and impact of the metallicity, [m/H], in the 
ATLAS 9 models requires special attention. Different metal- 
licities are achieved by scaling the abundances of all elements 
heavier than He by a single factor relative to the reference so- 
lar photosphere abundances ([m/H] = 0) given by Grevesse & 
Anders (1989). Clearly, a single parameter is insufficient to de- 
scribe the results of the complex physical processes which con- 
trol the variations in stellar abundance patterns. Nevertheless, it 
is equally clear that such a parameter can provide a useful diag- 
nostic if it can be determined robustly and with high precision 
from relatively low resolution data. 

The UV spectral region of all main sequence B stars 
(10,000 < T eff < 30,000 K) is extremely rich in Fe group 
lines (e.g., Walborn et al. 1995, Roundtree & Sonneborn 1993), 
which are obvious even at low resolution (Swings et al. 1973). 
The sensitivity of the ATLAS 9 models to changes in [m/H] 
is dominated by changes in Fe (except over narrow intervals 
where the presence of a few strong, isolated silicon and carbon 
lines becomes important, see Massa 1989). This assertion is 
justified in Figure [| where the solid curves show how of a -0.1 
dex change in [m/H] affects the UV SED's for three values of 
T e ff which span the range of the B stars. Below each curve, the 
shaded regions show the sum of the gf values (arbitrarily scaled 
and shifted) for the all spectral lines of Fe IV (top panel), Fe III 
(middle panel), and Fe II (bottom panel) within each ATLAS 9 
wavelength bin. Thus, the shaded regions crudely illustrate the 
relative strength and general wavelength distribution of the Fe 
opacity in the UV. In the figure, the Alog,F\ curves have been 
paired with Fe ions likely to be dominant, although absorption 
due to adjacent ionization stages will also be present. 

In each of the panels the opacity signature of the various Fe 
ions is strikingly similar to the signature of A [m/H]. It is clear 
that the value of [m/H] determined by fitting observed SED's 
with ATLAS 9 models is most heavily influenced by the abun- 
dance of Fe. Therefore we assume that the value of [m/H] de- 
rived from large scale UV absorption features in B stars is ef- 
fectively a measurement of [Fe/H]. 

Several comments regarding the robustness of an [Fe/H] 
measurement from fitting the UV SED's of the B stars can be 
made: 

1 . In addition to Fe itself, the other Fe-group elements (Cr, 
Mn, Fe, Co, Ni) also contribute to the UV opacity, albeit 
at a much reduced level. However, since all of these el- 
ements are thought to have similar nucleosynthetic ori- 
gins, their abundances are expected (and observed) to 
scale together naturally. Therefore the ATLAS 9 method 
of scaling all abundances by a single parameter is physi- 
cally reasonable for the Fe-group elements and does not 
confuse the relationship between [m/H] and [Fe/H] . 



2. The abundances of CNO and the light metals (e.g., Si) 
certainly have been observed to vary relative to Fe. For- 
tunately, however, the atmospheric opacity of main se- 
quence B stars is primarily due to hydrogen, helium and 
Fe. Consequently, even if an ATLAS 9 abundance mix 
results in too much or too little of CNO or the light met- 
als relative to the Fe group, the impact on the atmospheric 
structure and, hence, the shape of the low resolution en- 
ergy distribution is both minimal and localized in wave- 
length (see the discussion of Vega in §4) 

3. The low resolution UV spectral signature of Fe results 
from the combined opacity of literally thousands of lines. 
Therefore, random uncertainties in individual f-values or 
line centers tend to cancel out, as do non-LTE effects in 
specific lines. 

Thus, while the study of element-to-element abundance vari- 
ations generally requires high resolution stellar data and the 
computation of synthetic spectra, Fe abundances in main se- 
quence B stars can be determined readily from the gross struc- 
ture of their low resolution UV spectra. 

A second aspect of the ATLAS 9 abundances requires some 
clarification. The models use a reference solar Fe abundance 
which is 0.12 dex larger than the currently accepted value (see, 
e.g., Grevesse & Noels 1993). Since, as demonstrated above, 
Fe dominates the UV opacity signature from which [m/H] is 
determined, we adopt the following transformation: 



Fe 



H 



+ 0.12 



(5) 



where [m/H] is the metallicity parameter of the best-fitting AT- 
LAS 9 model and [Fe/H] is the actual logarithmic abundance of 
Fe relative to the true solar value. 

A final comment concerns the A [m/H] signature at the cool 
end of the B star range. The lower panel in Figure g (for 
T e ff = 10000 K) shows that, while the mid-UV correlation be- 
tween A log Fa and the Fe II opacity is unmistakable, the far- 
UV sensitivity to A [m/H] is larger than might readily be ex- 
plained by Fe alone. Part of this apparent discrepancy is due to 
changes in the atmospheric temperature structure which result 
from changing levels of opacity. The far-UV region in these 
relatively cool models is extremely sensitive to these changes. 
Thus, some of the rapid rise in A log F\ at short wavelengths 
is essentially a second order effect of changing the Fe opac- 
ity. The presence of strong opacity sources due to elements 
other than the Fe group also helps produce the far-UV sensitiv- 
ity to [m/H]. Numerous lines of C I are present at wavelengths 
shortward of ~4500 A, along with strong lines of C II, Si II, 
and other species. When this spectral region is included in the 
fitting process, the best-fit [m/H] value will represent a com- 
promise between the abundances of Fe and of these lighter el- 
ements. For the late B stars, an inability to achieve good fits 
between observed and theoretical SED's across the entire UV 
spectral region might serve to indicate a more complex abun- 
dance anomaly than can be represented by a single scale factor. 
In §4 we will demonstrate this effect in the spectrum of Vega. 

3. THE DATA 

The primary data used in our analysis are NEWSIPS low- 
dispersion IUE UV spectra (Boggess et al. 1978a, 1978b; 
Nichols & Linsky 1996), corrected and placed on the HST 



- FOS system (Bohlin 1996) using the algorithms developed 
by MF99. Data from both the long wavelength (LWR and 
LWP) and short wavelength (SWP) cameras are required and 
the wavelength regions covered are 1978-3000 A and 1 150— 
1978 A, respectively. For a few stars, UV-through-optical spec- 
trophotometry from the FOS is available. These data typically 
extend from 1 150 A into the optical spectral region. 

The UV/optical spectrophotometry (whether from IUE or 
FOS) are resampled to match the wavelength binning of the 
ATLAS 9 model atmosphere calculations in the wavelength re- 
gions of interest. Since the resolution of the IUE (~ 4.5 - 8 A) 
is only slightly smaller than the size of the ATLAS 9 UV bins 
(10 A), the characteristics of binned IUE data do not exactly 
match those of ATLAS 9, which may be thought of as result- 
ing from the binning of very high resolution data. In essence, 
adjacent IUE bins are not independent of each other, while ad- 
jacent ATLAS 9 wavelength points are independent. To better 
match the IUE and ATLAS 9 data, we smooth the ATLAS 9 
model fluxes in the IUE wavelength region by a Gaussian func- 
tion with a FWHM of 6 A. This slight smoothing simulates the 
effects of binning lower resolution data. No such smoothing 
is required when FOS data are used, since they have a much 
higher resolution than the IUE spectra. 

For FOS data, the weight assigned to each wavelength bin is 
derived from the statistical errors of the data in the bin. The 
weights for the IUE data are computed likewise, except that 
an additional, systematic, uncertainty of 3% is quadratically 
summed with the statistical errors in each bin to account for 
residual low-frequency uncertainties not removed by the MF99 
algorithms. 

We give zero weight to the wavelength region 1195 - 1235 
A due to the presence of interstellar Ly a absorption. Addi- 
tional points may be excluded in particular cases based on data 
quality or the presense of strong interstellar absorption lines 
(e.g., Fe II AA2344, 2383, 2600, Mg II AA2796, 2802) or stel- 
lar wind features (e.g., Si IV A1400, C IV A1550). Apparently 
spurious model features are present at 2055, 2325, and 2495 
A. These appear as isolated absorption dips in the models, but 
correspond to no observed stellar lines. Their strengths are tem- 
perature dependent, with the 2055 and 2495 A features present 
at model temperatures above ^30000 K and the 2325 A feature 
strongest at model temperatures below 25000 K (this feature is 
visible clearly in Figure 1, which also illustrates the metallic- 
ity dependence). We routinely give zero weight to these three 
wavelength bins. 

At present, we fit the satellite data and V magnitude data. 
Eventually, we will expand this to include Johnson UBV and 
Stromgren uvby photometry, once the calibration of the syn- 
thetic photometry is more secure (see §g). Our sources for 
the Johnson and Stromgren data are the catalogs of Mermilliod 
(1987) and Hauck & Mermilliod (1990), respectively, when 
available. 

4. RESULTS AND DISCUSSION 

In this section, we demonstrate our fitting process with sev- 
eral examples and discuss the consequences of each. The pro- 
gram stars are listed in Table [|. These stars were selected for 
three reasons. First, as a group, their physical parameters span 
the full range of the main sequence B stars. Second, funda- 
mental measurements exist for some of them, and these can be 
used to verify the accuracy of our results. Third, some of the 
stars have been the subject of fine analyses, so we can compare 



our derived stellar parameters to previous estimates, obtained 
by independent methods. 

BD+33° 2642 This is a high latitude post-AGB star and the 
exciting object of a planetary nebula. Its spectrum was exam- 
ined in detail by Napiwotzki, Heber, & Koppen (1994; hereafter 
NHK) who found T eff = 20200 K, logg = 2.9, and a peculiar 
chemical composition where C, N, O, Mg, and Si are subsolar 
by about 1 dex and Fe is subsolar by about 2 dex. This star is 
not typical of the type of object at which our research is aimed; 
however, it allows us to test the models in two respects. First, it 
is an extreme case of low metallicity. Second, the FOS data are 
of exceptional quality, allowing a critical assessment of how 
well the models can fit the details of an actual stellar energy 
distribution. 

We fit the FOS data for BD+33° 2642 using the R(V) = 3.1 
extinction curve from F99 and we solved for the 4 ATLAS 9 
model parameters, E(B-V), and the attenuation factor. The 
result is presented in Figure [| where we show the reddening- 
corrected stellar SED (small filled circles) overplotted with the 
best-fitting ATLAS 9 model (solid histogram-style curve). The 
fit is so good that the model and data are nearly indistinguish- 
able over most of the fitted wavelength range, with typical dis- 
crepancies less than 2% (consistent with the expected observa- 
tional errors) except in the regions indicated by crosses. These 
consist of the excluded ATLAS 9 bins discussed in §3 above 
and the regions containing the Ha and H/3 lines, which show 
evidence for emission in the FOS data. The model parameters 
are listed in columns 2-8 of Table g along with their la inter- 
nal fitting errors. The attenuation factor is represented by the 
computed angular diameter of the star 9 (in column 8, in units 
of milli-arcseconds), as noted in §2. 1 . 

Our temperature for BD+33 2642 is some 2000 K lower 
than that found by NHK. This is due entirely to the fact that 
NHK imposed a value of E(B - V) = 0.06, while our self- 
consistent value is significantly smaller (< 0.02). (The overcor- 
rection for extinction in the 2200 A region that results from as- 
suming E(B-V) = 0.06 is apparent in Fig. 3 of NHK). The other 
parameters derived by NHK will be adversely affected by the 
overestimated temperature, but are nonetheless generally in rea- 
sonable agreement with our results. Our value of [Fe/H] = -1 .5 
is consistent with the NHK Fe abundance ([Fe/H] = -2.0 ± 0.5 
from IUE data), but suggests that Fe may be more similar in 
abundance to the other elements than emphasized by NHK. We 
find a much lower value of v, than NHK, who adopted v t = 15 
km s" 1 based on line profile analysis. Such a large value is 
inconsistent with the observed UV opacity and might be an ar- 
tifact of their high T e jj. It would be interesting to repeat the 
NHK analysis using the lower temperature given in Table |[ 

HV 2274 in the LMC HV 2274 is a non-interacting eclips- 
ing binary system in the Large Magellanic Cloud which con- 
sists of two nearly identical early B stars. The properties of 
this system and its usefulness as a distance indicator for the 
LMC are discussed by Guinan et al. (1998). The FOS en- 
ergy distribution of HV 2274 (ranging from 1150 to 4800 A) 
was modeled using the program described here and the results 
can be seen in Figure 2 of Guinan et al. The fitting was per- 
formed by solving for the four model parameters, E(B-V), 
the attenuation factor, and the six additional parameters from 
FM which specify the shape of the UV extinction curve. The 
quality of the fit is comparable to that for BD+33 2642 dis- 
cussed above. The significance of HV 2274 for our purposes 
here is that the surface gravities of the stars in the system are 



extremely well determined from the properties of the light and 
radial velocity curves of the binary. This affords the oppor- 
tunity to test whether the values derived by our approach are 
reasonable. This is a critical test since it is not necessarily true 
that good fits to energy distributions produce the correct values 
of the stellar properties (e.g., see the discussion of HD 38666 
and HD 36512 below). The surface gravity of star A of the bi- 
nary system is known to be log g = 3 .54 ± 0.03 . When the fitting 
performed with gravity as a free parameter, a completely con- 
sistent value of \ogg = 3.48 ± 0.06 was determined. This result 
gives confidence in the values of the parameters and suggests 
that no large systematic effects are present. 

HD 38666, 36512, 55857, 31726, 3360, 61831, 90994, 
209952, 38899, 213398 These form a sequence of normal 
Pop I main sequence stars and are representative of the type of 
stars on which this study is mainly focused. This set of stars 
demonstrates the quality of fits possible across the entire range 
of temperatures typifying the B stars. In addition, HD 209952 
has a measured angular diameter (Code et al. 1976). Conse- 
quently, its attenuation factor (see eq. [Q]) can be compared to 
a direct measurement. 

The observational data for these stars are illustrative of that 
available for most of the objects to be studied (at least in the 
initial stages of the program) and consist of IUE UV spec- 
trophotometry (1150 - 3000 A) and Johnson and Stromgren 
optical photometry. The IUE spectra used here are listed in Ta- 
ble 3. The models which best reproduce the observed SED's 
are shown in Figure Q, overplotted on the dereddened and ar- 
bitrarily scaled observations. (The values of the scaling factors 
are given in the figure legend.) We show all the Johnson and 
Stromgren magnitudes (converted to flux units) for display pur- 
poses, but only the V magnitude was utilized in the fitting pro- 
cedure. Data points excluded from the fit (indicated by crosses) 
include the region of the strong Ly a interstellar line, the spuri- 
ous model points, and the C IV stellar wind lines in the hottest 
stars (see §3). As with BD+33 2642 above, we fit these lightly 
reddened objects by assuming the R(V) = 3.1 extinction curve 
from F99 and solving for the four model parameters, E(B - V), 
and the attenuation factor. 

While the data for these ten stars have lower S/N than for 
BD+33° 2642, the relative quality of the fits is the same. I.e., 
the discrepancies between the observed and model energy dis- 
tributions are all consistent with the expected uncertainties in 
the data (~3% for the IUE spectra and ~0.015 mag for the 
V magnitude). Note also that the measured value of the an- 
gular diameter 6 bserved f° r HD 209952 (last column of Table 
J2|) agrees with the value derived from the fitting procedure to 
within the expected errors. 

With the prominent exception of the two hottest stars, 
HD 38666 and HD 36512, the Fe abundances for the stars in 
Figure Qare within about ±0.2 dex of solar. The mean value for 
the eight stars with T eff < 30000 K is [Fe/H] = -0.07 ±0.15 
(s.d.). The dispersion is larger than expected from the mea- 
surement errors and may indicate a true range in the stellar Fe 
abundances. We will examine this further in the future with a 
larger data sample and look for corroborating evidence, such as 
a systematic dependence of abundance on spatial location. 

The results for HD 38666 and HD 36512 require some com- 
ment. The best-fitting models for these stars not only have 
very low (and very implausible) metal abundances, but also 
very high microturbulence velocities. While the combination 
of these two parameters yields models which reproduce the ob- 



served UV opacity features very well, nevertheless it seems 
likely that neither parameter actually represents a correspond- 
ing stellar property. 

We suggest that the large values of v, signify the presense of 
a non-thermal broadening which is not due to a random veloc- 
ity field (i.e., microturbulence) but rather to a velocity gradient 
in the photosphere — indicative of a stellar wind. The criti- 
cal difference, in terms of the effect on the line broadening, is 
that the continuity equation (for an atmosphere in whose den- 
sity decreases with altitude) requires that a velocity gradient be 
steepest in the outer layers of the atmosphere, while microtur- 
bulence is depth-independent in the ATLAS 9 models. Conse- 
quently, the large values of v, needed to fit the absorption from 
strongest lines (i.e., those formed in the outermost layers of 
the atmosphere) produce excess absorption by the less strong 
lines (formed deeper in the atmosphere). To compensate for 
this effect, the fitting program uses a lower metallicity model 
to weaken the latter, resulting in the unexpectedly low [m/H] 
values listed in Table @. If this explanation is correct then it 
signals a breakdown in one of the fundamental assumptions of 
the ATLAS 9 models, i.e., the applicability of hydrostatic equi- 
librium, and casts doubt on the correctness of the other model 
parameters despite the good fits to the observed SED's. We will 
examine this effect further with the larger sample, to delineate 
the region in the T e ff— logg plane where the values of v, begin 
to depart from the low values characteristic of most of the stars 
in Figure Q. 

HD 72660 and HD 172167 (Vega) In recent years it has 
been recognized that the "superficially normal" early A stars 
show a wide range of surface abundance peculiarities. HD 
72660 and Vega show distinctly different types of abundance 
anomalies and provide an excellent test of our ability to deter- 
mine Fe abundances. Furthermore, Vega, like HD 209952, has 
an angular diameter measurement listed by Code et al. (1976). 

Vega and HD 72660 have very similar temperatures and sur- 
face gravities. However, fine spectrum analyses of their op- 
tical Fe I and II lines have shown that Vega has a signifi- 
cantly sub-solar Fe abundance, [Fe/H] = -0.55 (Hill & Land- 
street 1993), while HD 72660 has a super-solar Fe abundance, 
[Fe/H] = +0.38, (Varenne 1999). Other elements show vary- 
ing degrees of departure from the solar standard. Carbon, for 
example, is nearly solar in Vega ([C/H] = -0.09) and strongly 
sub-solar in HD 72660 ([C/H] = -0.60). Only the coolest of 
the B stars have adequate optical Fe lines to allow accurate Fe 
abundances to be determined from the ground. Thus, verifica- 
tion of our approach to determining Fe abundances for Vega 
and HD72660 will add considerable confidence to our results 
for the hotter B stars, where we must rely exclusively on the 
UV. 

For Vega, we have both an excellent IUE data set (see Table 
||) and high precision optical ground-based spectrophotometry 
from Hayes (1985). The fit to Vega was performed as for the 
stars discussed above, by assuming the R(V) = 3.1 extinction 
curve from F99 and solving for the 4 ATLAS 9 model param- 
eters, E(B - V), and the attenuation factor. Both the UV and 
optical spectrophotometry (weighted using Hayes' errors) were 
incorporated in the fit and the ATLAS 9 models were smoothed 
by a 50 A gaussian in the optical region to mimic the resolution 
of the Hayes data. 

It turned out to be impossible to achieve a good fit to Vega 
simultaneously over the entire observed spectral region and, in 
the end, we gave the region 1200 < A < 1500 A zero weight. 



The results of this fit are shown as the top spectrum in Fig- 
ure H and the derived model parameters — the most notable 
of which is the low value of [m/H] = -0.73 — are listed in Ta- 
ble gj As is evident in the figure, the model which reproduces 
most of the observed spectrum to within the expected errors 
but significantly overestimates the output fluxes between 1200 
and 1500 A. The most likely reason for the poor fit in this re- 
gion is the complex abundance pattern in Vega. Because its 
Fe abundance is so low, the next largest opacity source in the 
far UV (C I and II) becomes important. The model shown in 
Figure g, which was optimized to fit the Fe-dominated mid-UV 
region, underestimates the C I and II opacity between 1200 and 
1500 A and thus overestimates the flux (see §2.6). The value of 
the Fe abundance for Vega derived from the fitting procedure, 
[Fe/H] = -0.61, agrees extremely well with the result from Hill 
& Landstreet (1993). A more detailed discussion of UV/optical 
energy distribution of Vega can be found in Castelli & Kurucz 
(1994). 

The fit to HD 72660 was performed in the same manner — 
but utilizing only the V magnitude in the optical region — and 
is shown as the lower spectrum in Figure || The fit is good 
throughout the entire range, although the overall quality of the 
observational data is lower than for Vega so that details in the 
residuals are not as apparent. The parameters of the fit are given 
in Table |2| notable is the large value of [m/H] = +0.25. Note 
that we did not need to disregard the shortest wavelengths in 
this case, as we did for Vega. The opacity is so high in this re- 
gion due to the elevated Fe abundance that the concurrent over- 
estimate of the C opacity has minimal impact. The derived Fe 
abundance for HD 72660, [Fe/H]= +0.37, is in excellent agree- 
ment with Varenne's result. 

HD 147933 (p Oph) This early B-main sequence star is 
well known for its peculiar UV interstellar extinction curve. 
We use this star to demonstrate the ability of our approach to 
extract information on the shape of the interstellar extinction 
curve simultaneously with the stellar properties. The observed 
SED of HD 147933 is shown in the top panel of Figure f| and 
consists of low-resolution IUE data (see Table |3|) and optical 
photometry. To fully constrain the models it was necessary to 
include not only the IUE data and the V magnitude in the fit, 
but also B — V, b—y, ml, cl, and u-b. Inclusion of the colors 
required a preliminary calibration of the synthetic photometry, 
which was derived from the differences between the observed 
and calculated photometry of the other program stars. We fit the 
HD 147933 data using the same six parameters as in the cases 
above plus we solve for the six additional parameters from FM 
which specify the shape of the UV extinction curve. Prior to 
the fit, the IUE data outside the range 1 195-1235 A (which are 
always given zero weight) were corrected for the effects of the 
strong damping wings of interstellar Ly a using a column den- 
sity of JV(HI) = 4.27 x 10 21 cm" 2 (Diplas & Savage 1994). 

The best-fitting model and the dereddened energy distribu- 
tion are shown in the middle panel of Figure £j The parame- 
ters of the model, along with E(B - V) and the angular diame- 
ter, are given in Table ||. The bottom panel shows the normal- 
ized extinction curve determined by the fitting program (thick 
solid line) overplotted on the actual normalized ratio of the 
best-fitting model to the observed (but Ly a-corrected) SED 
(filled circles). This result is compared with the "pair method" 
curve for HD 147933 published by FM (dotted line) and the 
F99 R = 3.1 curve (dash-dot line). The newly derived curve 
agrees with the published curve extremely well, lying within 
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the expected errors for a pair-method curve (Massa & Fitz- 
patrick 1986). This example demonstrates the feasibility of de- 
termining extinction parameters and stellar parameters simulta- 
neously. 

5. SUMMARY 

In this section, we summarize our results and expand upon 
several noteworthy aspects of the parameters listed in Table || 

1 . Our excellent fits to the entire sample of stars listed in Ta- 
ble |J which span the entire range of the main sequence B 
stars, are very encouraging. This is especially true in for 
the extremely detailed fit to the FOS spectrum of BD+33 
2642. 

2. Nevertheless, our results must be verified. Good fits do 
not necessarily mean that the derived physical parame- 
ters are correct (see item 4 below), and verification is 
essential. Therefore, it is gratifying that for those in- 
stances where we could compare the parameters derived 
from the models to directly measured quantities, the re- 
sults agreed. This includes the angular diameter measure- 
ments for Vega and HD 209952 and the surface gravity 
measurement for HV2274. 

3. Our Fe abundances for Vega and HD 72660 agree with 
previous measurements arrived at by completely inde- 
pendent means. This result gives us confidence that our 
metallicity determinations are actually faithful represen- 
tations of Fe abundances. It also justifies our interpre- 
tation of the general metallicity parameter in fits to the 
hotter B stars as Fe abundance. Furthermore the range in 
our observed Fe abundances is greater than the errors in 
the individual determinations. This suggests that study 
of a large sample of lightly reddened main sequence B 
stars might reveal patterns in the spatial distribution of 
abundances. 

4. Our result that v, ~ km s" 1 for all but the hottest main 
sequence B stars is comforting since it means that our 
findings do not have to rely on this poorly defined pa- 
rameter. It is also interesting that v, appears to assume 
non-zero values for the hottest stars in our sample. Per- 
haps this is verification of the long held suspicion that 
microturbulence is actually a surrogate for expansion in 
the early-type stars (e.g., Simonneau 1973; Massa, Shore, 
& Wynne 1992; Kudritzki 1992). This suggests that stars 
which require large values of v t may have atmospheres 
which are not in hydrostatic equilibrium. In these cases, 
the applicability of the Kurucz ATLAS 9 models is ques- 
tionable and the parameters derived from them are sus- 
pect. 

5. Our results for p Oph are very encouraging. They indi- 
cate that we may be able to extend our analysis in two 
directions. One is the determination of extinction curves 
for sight lines with far lower mean reddening than previ- 
ously possible, since our approach enables us to charac- 
terize the intrinsic flux of a star to far greater precision 



than the classic pair method approach, which typically 
relies on a coarse grid of standard stars. The other di- 
rection is the possibility of extending our investigations 
to include abundance determinations of more reddened 
stars, specifically those in the several open clusters ob- 
served by IUE. 

6. Our analysis also provides a powerful verification for the 
use of white dwarf models to calibrate UV instruments 
(Finley & Koester 1991, Bohlin 1996). Using a spec- 
trophotometric system calibrated by fitting white dwarf 
SEDs to model atmospheres, we obtained excellent fits to 
B star SEDs using an independent set of models whose 
input physics is very different. This amounts to a test 
of the internal consistency of the two sets of models and 
their abilities to reproduce observed SEDs. 

While the initial results presented here are extremely encour- 
aging, a few technical issues remain to be resolved. Leading 
among these are the following: 

1 . Verification of the calibration of the synthetic photometry 
computed from the ATLAS 9 models. 

2. The impact of binary companions on the parameters de- 
rived from the model fits. 

3. The impact of using the binned ATLAS 9 opacities to 
compute the theoretical SED's, instead of using high res- 
olution synthetic spectrum degraded to the same spectral 
resolution. 

4. The degree to which the precision of the results falls off 
with declining S/N in the observed SED's, as might be 
relevant for the application of this technique to extra- 
galactic stars. 

We will be examining these issues in future studies as we apply 
our results to the several hundred Milky Way stars for which 
the requisite data currently exist and begin to examine stars in 
Local Group galaxies. 

In summary, we have introduced and justified a procedure 
which employs UV spectrophotometry and optical photom- 
etry to simultaneously derive the intrinsic stellar parameters 
{T e ff,logg, [m/H],Vf} (where [m/H] is essentially [Fe/H]) of 
main sequence B stars and the impact of interstellar extinction 
on their spectra. We demonstrated that our approach produces 
excellent fits to the observational data and that reasonable phys- 
ical parameters are derived from the models as long as we allow 
the fitting procedure to access the full range of parameters and 
employ the FOS calibration for the UV data. This procedure 
represents a powerful astrophysical tool and mining the vast 
IUE archive of low dispersion B star data will be extremely re- 
warding. 

This program is supported by the Astrophysics Data Program 
grant NAG5-7113 to Villanova University and NASA contract 
NAG5-7372 to Raytheon ITSS. 
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Table 1 










Program Stars 






HD/BD 


Name 


Spectral 


V 


(B-V) 


(U-B) 


Number 




Type 


(mag) 


(mag) 


(mag) 


BD+33° 2642 


— 


B2IVp 


10.84 


-0.17 


-0.85 


HD 3360 


CCas 


B2IV 


3.66 


-0.20 


-0.85 


HD 31726 


— 


Bl V 


6.15 


-0.21 


-0.87 


HD 36512 


v Ori 


BOV 


4.62 


-0.26 


-1.07 


HD 38666 


/uCol 


09.5 V 


5.17 


-0.28 


-1.07 


HD 38899 


- 


B9IV 


4.90 


-0.07 


-0.17 


HD 55857 


- 


B0.5V 


6.11 


-0.24 


-0.99 


HD 61831 


- 


B2.5 V 


4.84 


-0.20 


-0.65 


HD 72660 


- 


Al V 


5.81 


0.00 


0.00 


HD 90994 


/3Sex 


B6V 


5.08 


-0.14 


-0.51 


HD 147933 


pOph 


B2IV 


4.59 


0.23 


-0.57 


HD 172167 


a Lyr 


AOVa 


0.03 


0.00 


-0.01 


HD 209952 


a Gru 


B7IV 


1.74 


-0.13 


-0.47 


HD 213398 


/3PsA 


A0V 


4.28 


0.01 


0.02 



Table 2 
Model Parameters 





Star Name 
(1) 


T eff 
(K) 

<2) 


log g 

(3) 


[m/Hl 1 
(4) 


[Fe/H] 1 

(5) 


n 

(tmr 1 ) 
(6) 


E(B-V) 

(mag) 
(7) 


e 2 

(mas) 

(8) 


a 3 

d/'.vcn'Ci/ 

(mas) 
(9) 



HD 38666 


31790 ±500 


3.85 ±0.08 


-0.74 ± 0.08 


-0.62 


14.7 ± 1.5 


0.021 ± 0.004 


0.10 ±0.01 


HD 36512 


31440 ±440 


3.88 ±0.21 


-0.82 ± 0.09 


-0.70 


1 1.0 ± 1.7 


11.031 ±0.004 


0.14 ±0.01 


HD 55857 


27740 ± 340 


3.83 ± 0.08 


-0.36 ± 0.07 


-0.24 


3.1 ± 1.0 


0.025 ± 0.004 


0.08 ±0.01 


HD147933 


24070 ± 460 


4.25 ±0.31 


-0.20 ±0.10 


-0.08 


0.0 ±0.0 


0.465 ±0.007 


0.42 ±0.01 


HD 31726 


23330 ±250 


4. 16 ±0.20 


-0.01 ±0.09 


+0.11 


1.3 ±0.8 


0.018 ± 0.005 


0.09 ±0.01 


HD 3360 


20900 ± 180 


3.44 ± 0.09 


-0.17 ±0.06 


-0.15 


2.7 ± 0.5 


0.013 ±0.004 


0.30 ±0.01 


BD+33°2642 


18170 ± 80 


2.45 ± 0.02 


-1.58 ±0.08 


-1.46 


6.3 ± 1.1 


0.012 ± 0.002 


0.01 ±0.00 


HD 61831 


17900 ± 140 


3.77 ±0.14 


-0.25 ± 0.09 


-0.13 


0.0 ±0.0 


0.017 ±0.004 


0.21 ±0.01 


HD 90994 


14420 ± 90 


3.64 ±0.11 


-0.16 ±0.05 


-0.04 


0.0 ±0.0 


0.006 ± 0.004 


0.22 ±0.01 


HD209952 


13440 ± 50 


3.76±0.11 


-0.35 ± 0.05 


-0.23 


0.6 ±0.7 


0.000 ± 0.000 


1.07 ±0.01 1.02 ±0.07 


HD 38899 


1 11)4!) - 20 


3.80 ±0.06 


-0.20 ± 0.03 


-0.08 


1.1 ±0.5 


0.000 ± 0.000 


0.30 ±0.01 


HD2 13398 


9660 ± 20 


4.06 ±0.04 


+O.08 ± 0.02 


+O.20 


0.9 ±0.4 


0.000 ± 0.000 


0.45 ±0.01 


HD 172 167 


9600 ± 10 


3.84 ±0.04 


-0.73 ± 0.06 


-0.61 


1.9 ±0.4 


0.000 ± 0.000 


3.26 ±0.01 3.24 ±0.07 


HD 72660 


9510 ± 20 


3.95 = 0.05 


+0.25 ± 0.03 


-K).37 


2.8 ±0.3 


0.000 ± 0.000 


0.11 ±0.01 



The quantity [m/Hl in column 4 is ihc "melallieily" of the best-filling ATLAS 9 model. [Fe/H] in column 5 represents the relative abundance of Fc with respect 
to the currently accepled solar value, compuled from [in/H] via equation 5. See the discussion in §2.6 Nine dial [lie results for the- sLars lit) 38666 and I ID 36512 are 
unlikely to icpresenl llie true mendicities of these stars, but rather rcsull from complicating effects due to the onset of strong stellar winds. See the discussion in S4- 

In columnS, the scale factors (Rid) determined from the fits have been converted into angular diameters (in units of milli-arcseconds) as noted in §2.1. 

The entries in column 9 are directly measured stellar angular diameters from Code el al. 1976. 



11 



Table 3 
IUE Spectra Used for Program Stars 



HD 


SWP image 


LWP image 


LWR image 


Number 


numbers 


numbers 


numbers 


HD 3360 


4316,26364,26365,26366, 


2931,4866,5041,5042, 


3812 




26376,26491,26493,26535, 


6503,6505,6507,6570, 






26647,26873,26874,26875, 


6692,7528,7529,7530, 






32869,37984 


12616,17121 




HD 31726 


8165 




7098 


HD 36512 


8164,11155 




7054,7097,9787 


HD 38666 


14340 




10954 


HD 38899 


16639 




12875 


HD 55857 


14339,23117 


3449 


10953 


HD 61831 


14309 




10940 


HD 72660 


44438 


22854 




HD 90994 


7927,9219,15791 




6905,7975,12162 


HD 147933 


6587,6588,6589 




5639,5640 


HD 172167 


27024,29864,29867,30548, 


7010,7038,7888,7889, 






30549,32906,32907 


10347,10348 




HD 209952 


26995 


7009 




HD 213398 


42704 


21480 
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0.02 

0.01 
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0.01 

0.0 

-0.01 

0.01 
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-0.01 

0.01 
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-0.01 

-0.01 



&£ -0.02 
o 



-0.03 



T = 20,000 K 

eff 

log g = 4.0 
~m/H] = 

2.0 km/s 





AT = +100 K 



Alog g = +0.1 dex 




AE(B-V) = +0.005 



1500 



2000 2500 3000 3500 
Wavelength (A) 



4000 



4500 



FIG. 1. — Representative spectral signatures of the four model parameters {T e ff.\ogg,[m/H],Vt} and E(B-V). The top panel shows the normalized flux for a 
model with T t ft = 20000 K, logg = 4, [m/H] = 0, and v t = 2 km s~' . Each of the lower panels shows the logarithmic change in the surface flux of the model when 
one of the parameters is changed by the amount A indicated in the panel. The details of the spectral signatures, except that for E(B — V), are dependent on the basic 
parameters of the models, particularly on T e tf. 
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FIG. 2. — The spectral signature, Alog/^, due to a change of A[m/H] = —0.1 dex at three different values of T e rf (solid curve in each panel), compared with 
a measure of the strength of the absorption due to various ions of Fe (shaded regions). In particular, the shaded regions show the total values of gf for all the 
Fe transitions within each ATLAS 9 wavelength bin, arbitrarily scaled and shifted for display. The figure demonstrates that, in the UV region and within the 
temperature range of the B stars, the overall signature of changes in the "metallicity" of an ATLAS 9 model is mainly due to one particular element — Fe. Other 
elements contribute within much more restricted wavelength regions, such as Mg II A2800 (in the T e et = 10000 K case) or Si IV A1400 and C IV A1550 (both in the 
T e ff — 30000 K case). The strong features in AlogF^ at 2325 A and 2495 A in the 20000 K and 30000 K panels, respectively, are spurious model features which 
correspond to no observed stellar lines (see §3). 
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FIG. 3 .— The dereddened FOS energy distribution of BD+33° 2642 from 1 150 A to 8000 A is shown (solid circles) along with the best-fitting ATLAS 9 model 
(solid histogram line). Crosses indicate data points excluded from the fit, due to contamination from interstellar gas-phase absorption lines, spurious model points, 
and emission in the stellar Ha and H/3 lines. The inset shows an expanded view of the region surrounding the Balmer jump. The parameters of the model are given 
in Table 0; the R(V) = 3.1 extinction curve from F99 was used. Over the wavelength range shown, the deviations between the model and the data are typically less 
than 2%. 
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FIG. 4. — Sample fits for lightly reddened stars. Dereddened energy distributions for 10 main sequence stars ranging from spectral type 09.5 to A0 V are shown 
(symbols) along with the best-fitting ATLAS 9 models (solid histogram lines). Crosses indicate data points excluded from the fits. The parameters of the models 
are given in TableH; the R(V) = 3.1 extinction curve from F99 was used. The Johnson and Stromgren photometry data have been converted to flux units for display 
only. The fitting procedure utilizes only the UV fluxes and the V magnitude. The U band point appears high in stars with large Balmer jumps because its broad band 
straddles the jump. The logarithmic fluxes have been arbitrarily displaced vertically for display purposes. The offsets applied are, from top to bottom, 14.22, 13.59, 
13.80, 13.41, 12.03, 12.08, 11.78, 10.06, 10.90, and 10.26 
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FIG. 5. — Best-fitting ATLAS 9 models for HD 172167 (Vega; top) and HD 72660 (bottom). The symbols are the same as in Fig. Wand the same reddening model 
was used. The fit for Vega is relatively poor between 1200 < A < 1500 A, possibly due to carbon opacity unaccounted for by the model (sealext). Both stars have 
similar temperatures and surface gravities, and their energy distributions differ primarily because of their different Fe abundances (see Table 0). 
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FIG. 6. — Best-fitting ATLAS 9 model and UV/optical extinction curve for HD 147933 (p Oph). The top panel shows the observed IUE plus optical energy 
distribution (symbols same asJn Fig. 4). The middle panel shows the dereddened energy distribution along with the best-fitting ATLAS 9 model. The parameters of 
the model are given in Table pj. The dereddened energy distribution was also corrected for the effects of interstellar Ly a absorption, assuming a column density of 
N(W) = 4.27 X 10 2 ' cm~ 2 from Diplas & Savage 1994. The lower panel shows the normalized UV extinction curve determined by the fitting procedure (thick solid 
line) overplotted on the actual normalized ratio of the best-fitting model to the observed (but Ly a-corrected) SED (filled circles). For comparison, the normalized 
extinction curve for HD 147933 derived by FM using the traditional pair method (dotted line) and the average Galactic curve for R(V) = 3.1 from Fitzpatrick 1999 
(dash-dot line) are also shown. 



